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INTRODUCTION 


Fatigue  fractures  occur  at  points  of  stress  concentration, 
which  are  normally  regions  of  complex  stress  distribution  with 
large  gradients.  The  stress  field  is  geometry  dependent  and 
will  only  be  reproduced  in  a  laboratory  sample  if  geometric 
similitude  is  preserved.  Matching  sample  to  part  must  also 
include  similitude  in  the  two  major  phases  of  fatigue;  initia¬ 
tion  and  propagation  and  not  merely  a  matching  of  stress  con¬ 
centration  factors.  In  fact  if  the  fatigue  process  were  well 
understood,  it  would  seem  that  a  knowledge  of  the  local 
conditions,  coupled  with  material  property  information,  would 
be  all  that  is  required  to  predict  fatigue  life. 

This  has  not  been  possible  so  far,  but  this  study  has  been 
undertaken  to  see  how  much  success  can  be  achieved  toward  this 
end  via  a  detailed  study  of  conditions  surrounding  the  notch 
tip.  To  have  any  chance  of  success  whatsoever,  the  local  stress 
and  strain  conditions  at  the  notch  must  be  determined  accu¬ 
rately.  One  of  the  means  available  for  accomplishing  this 
is  experiment,  and  this  approach  was  used  successfully. 

SIMILITUDE 

An  evaluation  of  how  well  one  notch  specimen  models 
another  can  be  studied  via  similitude  theory.  The  strains  in 
the  vicinity  of  a  notch  depend  upon  the  following  variables: 

£(.  :  +i(  ?,  W,  *1/9,  (1) 
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where 


X*%  —  Strain  components 
P  =  Applied  load 
W  =  Sheet  width 
t  =  Sheet  thickness 
^  -  Radius  of  curvature  of  notch 
d  =  Depth  of  notch 
E.  =  Tensile  Modulus 
0  =  Poisson's  Ratio 

By  dimensional  analysis,  the  functional  relationship  can  be 
written  in  terms  of  groupings  of  variables  that  are  dimen¬ 
sionless  . 


&iz  -K( f  >  y1  4’  4 >°) 


(2) 


where 
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From  similitude  theory,  if  in  two  different  geometries 


Then 


GEOMETRIC  SIMILARITY 


If  the  material  of  the  part  and  the  model  are  the  same, 
the  same  strains  are  produced  only  if  they  are  geometrically 
similar.  It  is  not  sufficient  to  match  only  the  stress  at 
the  notch  tip,  but  complete  similarity  must  be  maintained. 

If  the  complete  strain  fields  are  not  the  same,  the  plastic 
zones  will  be  different,  and  the  neighborhood  where  fatigue 
action  is  taking  place  has  different  conditions.  Thus, 
similarity  must  be  preserved  when  employing  experimental 
fatigue  results  from  one  specimen  to  predict  fatigue  failure 
in  another.  This  is  not  generally  done  in  current  practice. 
Since  it  is  not  convenient,  nor  even  practical,  to  run  tests 
on  a  geometrically  similar  specimen  each  time  a  part  is  an¬ 
alyzed  for  fatigue,  a  different  approach  is  attempted  to 
develop  a  more  general  theory  where  knowledge  of  local 
conditions  at  the  notch  tip  is  sufficient,  along  with 
material  properties,  to  predict  fatigue  in  any  specimen. 

SIMILITUDE  IN  UNIAXIAL  TENSION 

Because  much  of  the  fatigue  data  available  has  been 
collected  on  uniaxial  tension  specimens  and  because  some¬ 
thing  resembling  uniaxial  tension  prevails  at  notch  tips, 
it  is  reasonable  to  pursue  this  similarity,  which  is  shown 
in  Figure  1. 
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Figure  1.  Uniaxial  Model  of  Stress  Condition  at 
a  Notch  Tip. 

For  an  elastic  material 

ii  -  'k(S'*,cry>crti  (3) 

If  the  same  material  is  used  for  the  model  and  the  specimen 
the  strain  field  will  be  modeled  to  the  degree  that  the 
stress  field  is  modeled.  For  instance,  in  the  region  of 


t ne  notch  tip  in  plane  stress, 

» -L  [5  -  0<ry]  (4, 

To  keep  the  departure  from  a  uniaxial  stress  state  less  than 
3%,  which  would  be  in  the  range  of  experimental  measurement 
error,  the  transverse  stress  must  be  small. 

(J^  <  0.05  (5) 

For  a  typical  value  of  Poisson's  Ratio  for  aluminum,  "0  ■  0.33 
and 

(Tx  ^  O.  °  '  ^*y 

In  the  y  direction. 

To  maintain  the  same  3%  limit, 

■OC;  a  o.05  Gy 

or 

(j;  <  0.09  <jy 

Also  in  the  Z  direction, 

C%  i  o-ozO) 


(6) 


(7) 


(8) 


(9) 
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The  most  demanding  of  these  is  the  first.  In  order  to 
model  the  behavior  at  the  notch  tip  as  a  uniaxial  specimen, 
only  regions  may  be  considered  where  the  radial  stress  is  one 
per  cent,  or  less,  of  the  hoop  stress. 

For  a  circular  hole  in  an  infinite  sheet,  the  location 
where  -  0.01  is  at  Y"  =  1.0101  .  For  a  two 

inch  diameter  hole,  the  distance  from  the  edge  of  the  hole 
where  the  error  exceeds  the  stated  allowable  is  0.01  in. 

This  precludes  the  use  of  a  strain  gage  to  make  local  measure¬ 
ments.  If  specimen  sizes  were  enlarged  ten  times  until  the 
uniaxial  tensile  region  were  large  enough  to  accommodate  a 
1/8"  strain  gage,  the  hole  would  be  twenty  five  inches  in 
diameter  and  the  specimen  would  be  four  to  six  feet  wide  and 
ten  to  twelve  feet  long.  To  circumvent  this  difficulty, 
photoelasticity  was  used  as  the  experimental  method  for 
measuring  local  strain  conditions  at  the  notch.  There  are 
similar  si2e  problems,  but  they  are  not  quite  as  severe  as 
for  strain  gages,  plus  the  bonus  of  obtaining  the  entire 
strain  field  at  the  notch  tip  made  photoelasticity  attractive. 

PHOTOELASTIC  MATERIAL  PROPERTIES 


SELECTION 

7075-T6  aluminum  was  selected  for  the  study  of  notch 
fatigue  since  it  is  a  popular  alloy  of  the  aircraft  industry 
and  one  for  which  there  is  a  sizeable  fatigue  data  base. 
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To  coat  7075-T6  with  a  birefringent  material  that  may  be 
loaded  into  the  plastic  zone,  Photolas tic* s  PS-1  material  was 
selected.  It  has  a  maximum  strain  capability  of  10%  and  is  a 
good  all  around  photoelastic  material. 

Mil  Handbook  5  gives  typical  data  for  7075-T6  sheets 
of  0.080  in  thickness.  To  extend  into  the  yield  range,  strains 
of  10,000  ps  must  be  measured  with  excursions  up  to  20,000  ps 
probably  anticipated.  From  equation  (A-l)  of  Appendix  A,  it 
is  readily  seen  that 


If  something  like  five  fringes  are  desired  at  the  10,000  ps 
level,  the  F^  value  needed  is  0.002.  This  would  allow  10 
fringes  at  the  maximum  strain  level  expected. 

From  equation  (A-4) ,  the  material  fringe  value  is  given 


K 


(id 


For  PS-1,  the  company  supplied  a  value  of  K  =  0.15,  and  for 

white  light  the  wave  length  is  ^  =  22.7  x  10~6  in;  thus 

-4 

the  material  fringe  value  is  1.5  x  10  .  Using  equation  (A-l) , 

the  thickness  of  the  coating  desired  is 

V\t  T  _  1,5  *  _  -  0.0 56 

ZFt  U)(-ooi) 
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The  coating  material  PS-lC  was  purchased  in  sheets  nominally 
.040  in  thick.  The  bonding  agent  used  was  that  suggested  by 
Photolastic.  At  first  PC-1  was  used,  but  for  higher 
elongations  it  was  necessary  to  go  to  PC- 8,  even  though  the 
first  cement  was  supposed  to  have  the  elongation  capability 
required. 

CHARACTERIZATION  OF  PS-lC  PHOTOELASTIC  MATERIAL 

A  careful  characterization  of  PS-lC  was  conducted  by 
Stenstrom  (1)  using  two  uniaxial  tensile  test  specimens 
prepared  from  a  single  sheet  of  the  photoelastic  material. 

The  first  specimen  was  1.0  inch  wide  and  was  used  for  thick¬ 
ness  measurement  and  longitudinal  strain  data.  The  second 
2.0  inch  wide  specimen  was  used  for  transverse  strain  data 
and  for  measuring  the  photoelastic  coating  stress  fringe 
value,  Fg.  . 

Thickness  Determination 

One  longitudinal  edge  of  the  1.0  inch  wide  specimen  was 
lightly  sanded  to  create  a  smooth,  uniform  cross-section  of 
the  photoelastic  material.  It  was  then  mounted  such  that  a 
traveling  microscope  could  be  used  to  measure  thickness.  With 
the  contrast  in  appearance  between  the  reflective  coating  and 
the  base  material,  it  was  possible  to  accurately  determine 
the  thickness  of  each.  For  the  specimen  used,  the  base 
thickness  was  measured  to  be  0.0382  inches  and  the  reflective 
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coating  to  be  0.0051  inches.  These  numbers  were  obtained  by 
measuring  eight  separate  points  along  the  edge  and  averaging 
the  results.  The  total  average  thickness  of  0.0433  inches 
fell  within  the  range  of  the  manufacturer's  measurement  of 
0.044+0.002  inches.  It  was  therefore  concluded  that  the 
labeled  thickness  included  both  the  base  and  the  reflective 
coating  after  two  similar  specimens  of  different  marked 
thicknesses  were  checked.  For  calculation  purposes,  the 
thickness  of  the  coating  was  taken  to  be  0.0382  inches  since 
the  path  of  polarized  light  does  not  penetrate  the  reflective 
coating  of  the  photoelastic  material.  Inherent  in  using 
this  thickness  is  the  assumption  that  the  reflective  coating 
does  not  contribute  to  the  reinforcement  of  the  base  material. 
Young  1 s  Modulus 

Two  lines  were  lightly  etched  on  the  surface  of  the  1.0 
inch  wide  specimen  in  the  horizontal  direction  in  order  to 
accommodate  a  Tuckerman  Optical  Strain  Gage.  The  specimen 
was  then  loaded  in  tension.  The  readings  were  taken  with  a 
autocol lima ter  and  recorded.  Using  the  American  Instrument 
Company  calibration  factors,  the  strain  was  calculated  with 
the  following  formula  taken  from  Ref.  2. 


S  =  P'A'L,R- 
\oooe 


(12) 
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where  S 


=  Strain  in  microinches  per  inch 
F  *  Lozenge  calibration  factor,  =  1.002 
A  =  Auto-collimater  calibration  factor,  =  1.004 
L  =  Lozenge  length,  2.0  inches 
R  -  Net  reading 
e  =  Gage  length,  1.0  inch 

With  the  measured  strain,  the  measured  load  at  each  point 
and  the  known  cross-sectional  area,  Young's  Modulus,  E,  was 
determined  from  the  slope  of  the  stress-strain  curve  shown  in 
Figure  1,  using  a  linear  regression  analysis. 

<Tz  3S7,Ol“i.'S  6,  -  4  *“>7  <i3) 

R.5  0.9968 

where  R  *  correlation  coefficient. 

Dropping  the  non-zero  intercept.  Young’s  modulus  was  taken 
to  be 

Et  =  pSL  (14) 

Photolastic  Inc.  information  lists  this  modulus  to  be  nomin¬ 
ally  360,000  psi.  Using  thin  specimens  for  this  measure¬ 
ment  requires  very  accurate  thickness  measurements  because 
E  is  inversely  proportional  to  thickness.  Small  errors  in 
thickness  can  turn  out  to  produce  large  percentage  errors  in 
area  and  modulus. 
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PS-1C  CHARACTERIZATION 
STRESS-STRAIN  CURVE 


Figure  1 


Poisson's  Ratio 


Using  the  2.0  inch  wide  PS-1C  specimen,  two  lines  were 
lightly  etched  into  the  surface  in  the  longitudinal  direction 
one  inch  apart.  With  the  Tuckerman  Optical  Strain  Gage 
mounted  in  the  transverse  direction  as  shown  in  Figure  2, 
the  specimen  was  loaded  in  tension.  The  data  points  from 
the  test  were  used  in  a  linear  regression  analysis  with  the 
following  results: 

(y  -  £2.  -  I.  ooio&  (is) 

R. C  0.99991 

Poisson's  ratio  was  then  determined  by  dropping  the  non-zero 
intercept  and  dividing  the  negative  of  the  slope  by  E  as 
indicated  in  Figure  1.  The  measured  Poisson's  ratio  is 

0  =  0 .  3793 

v 

which  compares  with  a  nominal  value  of  0.38  supplied  by 
Photolastic  Inc.  These  tests  were  performed  in  accordance 
with  ASTM  standards  for  uniaxial  properties  (3) . 

Coating  Fringe  Value  in  Terms  of  Stress,  F0 

Prior  to  etching  lines  on  the  2.0  inch  wide  specimen,  it 
was  loaded  in  tension  in  the  elastic  region,  and  the  fringe 
number,  N,  was  determined  using  a  reflection  polariscope. 

With  several  load  points  and  known  dimensions,  the  stress 
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re  2.  TUCK  HUMAN  OPTICAL  STRAIN  GAGE  SLT-IJI 


was  calculated  and  plotted  against  the  fringe  number.  Using 
the  data  from  Table  I,  a  linear  regression  was  made  with  the 
data  which  yielded 


CTz  A-Sio.UU  N  +  1.9114  U6) 

Dropping  the  non-zero  intercept,  the  stress  fringe  value  was 
taken  to  be 

(17) 

Coating  Fringe  Value  in  Terms  of  Strain,  Fe 

With  the  tensile  modulus  and  the  Poisson's  Ratio  of 
the  coating  material  already  determined,  the  strain  coating 
fringe  value  can  be  calculated  readily.  From  Appendix  A, 

Ft  =  t  F<j*  -  0"U  l*/c*/-f 

E.C. 

CHARACTERIZATION  OF  7075-T6 
The  elastic  and  plastic  material  properties  of  the 
aluminum  panels  were  established  by  KAISER  (4)  using  tensile 
tests  of  uniaxial  specimens  made  from  the  same  mill  run  as 
the  sheets  to  be  tested.  The  specimens  were  manufactured 
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PS -1C  I'  ,  DETERMINATION 


u 

Compensator 

psi 

Reading 

N 

0.0 

0.0 

0.0 

202.37 

23.0 

0.489 

408.42 

39.0 

0.830 

610.26 

63.0 

1.340 

785.26 

79.5 

1.691 

990.53 

102.0 

2.170 

1189.21 

123.0 

2.617 

LINEAR  REGRESSION : 


and  tested  according  to  ASTM  standards  E8  and  B557.  (3)  The 
specimens  were  cut  transverse  to  the  rolling  direction,  and 
the  tests  were  performed  in  a  100  Kip  MTS  machine.  The 
duration  of  the  tests  was  three  minutes.  MICRO-MEASUREMENTS, 
EA-13-125AD-120 ,  precision  strain  gages  with  a  temperature 
compensated  bridge  circuit  were  used  on  all  specimens.  Trans' 
verse  gage  sensitivity  errors  were  corrected  according  to  the 
manufacturer's  recommendations  (5).  Critical  cross-section 
measurements  were  made  with  a  micrometer. 

YOUNG'S  MODULUS  (E) 

Tests  were  conducted  using  the  specimen  shown  in 
Fig.  3  on  the  MTS  test  machine  with  10,000  lb  INSTRON  grips, 
which  gripped  the  specimen  evenly.  The  results  of  testing 
three  specimens  are  shown  in  Tables  II  to  IV  where  all  strain 
values  have  been  corrected  for  transverse  sensitivity. 

Linear  regression  of  the  first  eight  points  in  the  elastic 
range  of  all  three  specimens  gave 

G*  =  +■  l*Ol»7  X  ID1  t  (18) 

2 

with  a  correlation  coefficient  of  R  *  0.9997.  Ignoring 
the  non- zero  intercept  as  specified  by  ASTM 

£  r  I  0.  I  2.  X  lOfc  psi.  (19) 
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TABLE  II.  MTS  SPECIMEN  A  TEST  RESULTS 


Cross-section  =  0.03975  in2 
Load,  lbs. _ _ Strain,  e,10~6  in/in 


256 

*-*«•****  f  C-  f  Jm  \j  J-H/  Xll 

615 

503 

120  4 

750 

1801 

1005 

2423 

1255 

3042 

1505 

3665 

1778 

4352 

2003 

1 

4925 

2252 

5587 

2508 

6355 

2755 

7365 

2905 

8230 

2984 

9045 

3037 

10150 
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TABLE  III.  MTS  SPECIMEN  B  TEST  RESULTS 


2 

Cross-section  =  0.03975  in 
Load,  lbs. _ Strain , 


500 

1250 

750 

1860 

1000 

2450 

1250 

3080 

1500 

3700 

1750 

4320 

2000 

4950 

2250 

5630 

2500 

6380 

2750 

7340 

2900 

8200 

3000 

9250 

3050 

10200 

3100 

11550 

3125 

12  800 

,  10  f  in/in 


1 
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TABLE  IV.  MTS  SPECIMEN  C  TEST  RESULTS 


Cross - 

section 

=  0.03975  in 

Load, 

lbs . 

Strain 

272 

663 

503 

1220 

762 

1850 

100  8 

2  451 

12  31 

3005 

1506 

3690 

1755 

4318 

2000 

49  40 

2255 

5607 

2503 

6325 

2750 

72  30 

2900 

8120 

3000 

9200 

3060 

10500 

3100 

11500 

3125 

12500 

,  10  ^  in/in 
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POISSON'S  RATIO  (v) 

Tests  were  conducted  using  the  specimen  shown  in 
Fig.  4  on  the  RIEHLE  test  machine  with  10,000  lb  RIEHLE 
grips.  The  results  were  tabulated  in  Table  V.  The  trans¬ 
verse  strain  measurements  were  corrected  by  the  manufacturer's 
data  on  transverse  sensitivity.  Both  strains  were  read  at 
the  same  load  and  plotted  versus  one  another.  A  linear  re¬ 
gression  was  made  on  the  data  resulting  in: 

£  r  0.0000  —  O.bZSt*  t. ,  (20) 

with  a  correlation  coefficient  of  R  =  0.99998.  Disregarding 
the  constant  factor  as  insignificant  to  four  places 

■0  =  O.iZSk  (2i) 

YIELD  STRESS  AND  STRAIN  HARDENING  MODULUS 

The  values  of  yield  stress  and  strain  hardening  modulus 
were  determined  graphically  using  the  data  from  Tables  III 
and  IV.  Plastic  region  data  in  Table  II  is  not  reliable 
because  of  excessive  creep  encountered  during  that  test. 

0.2%  offset  yield  stress,  =  76,000  psi 
strain  hardening  modulus,  E.  =  566,000  psi 
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FIGL’PE  4 
GAGE  SPECIMEN 


AXIAL  STRAIN 
GAGE 


TRANSVERSE 
STRAIN  GAGE 


thicknctt  =  0.0800'* 


- 1.5  -*-j 
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TABLE  V.  REIHLE  SPECIMEN  TEST  RESULTS 


2 

Cross-section  =  0.12  in 


LOAD 

lbs. 

STRAIN, e  x 

10  ^in/in 

STRAIN,e2 
10  ® in/in 

500 

345 

-115 

1000 

720 

-240 

1500 

1114 

-364 

2000 

1501 

-492 

2500 

1895 

-614 

3000 

2291 

-745 

3500 

2705 

-871 

4000 

3095 

-998 

4500 

3515 

-1126 

5000 

3912 

-1255 

5500 

4332 

-1384 

6000 

4750 

-1520 

6500 

5188 

-1647 

7000 

5595 

-1784 

7500 

6075 

-1928 

8000 

6649 

-2108 

8500 

7385 

-2375 

9000 

8663 

-2888 

9500 

12245 

-4435 
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The  graphical  fit  of  these  values  to  the  test  data  can  be 
seen  in  Fig .  5 . 

RAMBERG-OSGOOD  COEFFICIENTS 

The  Ramberg-Osgood  equation  for  elastic-plastic 
stress-strain  characterization  is  given  by: 


where : 


£  =  strain 
(T  =  stress 
£  =  Young's  modulus. 

The  6  and  n  coefficients  were  determined  graphically  from 
the  data  of  Table  IV,  by  the  method  given  by  Ramberg  and 
Osgood  (6) . 


.n-i 

z3  1  7  [  ~k] 


(23) 


Y\  =  ) 


(,7A) 


(24) 


Lines  corresponding  to  slopes  of  0 . 7E  and  0 . 85E  were 
laid  out  on  the  X-Y  recorder  plots  of  load  versus  longitudinal 
strain.  Load  values  where  the  rays  intersected  the  load- 
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strain  curve  were 


P  =  3060  lbs  P  =  2950  lbs 

»1£  *• 


The  cross-sectional  area  from  Fig.  3  is  0.03975  in  .  This 
gave 

77,484  psi.  Gm  :  14,7-14  p*. 


The  Ramberg-Osgood  Parameters  are  given  by 

fa  r  1.475  x  104* 

rv  -  zi.se> 

The  stress-strain  law  for  uniaxial  tension  throughout  the 
entire  range  is 

-8^.  „  *  -  io3  —  XI. Sd 

9.86HU0  V  +.9949KiO  <T  (25) 

This  is  shown  plotted  in  Fig.  5  along  with  data  points  from 
two  of  the  specimens. 


REINFORCEMENT 

When  specimens  are  made  from  thin  sheets  of  material 
to  duplicate  wing  skins,  the  coating  material  can  provide 
some  reinforcement  to  the  specimen  and  influence  the  stress 
values.  This  is  usually  small,  but  since  it  is  calculable 


by  classical  developments,  allowance  can  be  readily  made 
for  it. 


Consider  an  infinitesimal  element  in  a  biaxial  stress 
field  oriented  along  principal  stress  directions  as  illus¬ 
trated  in  Fig.  6.  If  the  total  loads  applied  to  the 
element  are  P^  and  P£,  they  can  be  equated  to  the  loads  in 
the  coating  and  the  specimen, 

P,  =  P,4  *  P.c  Pj  --  P.»  ♦  P«.  <; 

If  there  were  no  coating,  the  stress  produced  by  and  P2 
would  be 

P,  =  (Tu  k*  1  ( 21 

The  stresses  in  the  component  parts  of  the  coated  specimen 
are  given  by 

Pss-  (28 

P«-s  <r“W‘<ix'  <29 

Equating  the  loads  in  the  coated  and  the  uncoated  elements, 

cr,K  =  +  crickc.^u*.  (3C 
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FIGURE  6.  Reinforcement  Model 
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Simplifying  terms. 
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The  quantities  in  the  bracket  have  all  been  measured;  there¬ 
fore  all  photoelastic  measurements  must  be  multiplied  by 
this  factor  to  obtain  uncoated  specimen  response 

-  1.0  H,1  ftM  (39) 

POISSON'S  RATIO  MISMATCH 

As  mentioned  previously  the  notch  tip  may  be  considered, 
under  certain  restrictions,  as  a  uniaxial  specimen.  As 
the  tip  is  strained,  the  lateral  contraction  is  different  in 
the  specimen  and  in  the  coating.  At  the  interface  of  the 
coating  and  specimen, 

Uc  =  =  ’  ^  =  '  ^  t'C  (40) 

at  the  free  surface  of  the  coating 


Thus  there  is  a  transverse  strain  gradient  through  the  thick¬ 
ness  of  the  coating.  The  birefringence  will  exhibit  the 
average  strain  readings,  which  will  lie  somewhere  between  the 
values  at  the  interface  and  the  values  at  the  surface.  At 
the  interface. 


(42) 


At  the  free  surface, 


tic'  ^ 


tc- 


-  €..c  C14^)  ' 


(43) 


The  actual  fringe  value,  N,  will  be  an  average  bounded  by 
these  two, 

6tc.  ( H-0»)  <.  N  4,  (ltOc) 

Ft 


(44) 


Experiments  with  uniaxial  tensile  specimens  have  shown  (7) 
that  the  strain  variation  through  the  thickness  are  only 
present  near  free  boundaries  in  the  plane  of  the  specimen, 
such  as  the  notch  surface.  This  arises  because  of  the  edge 
effect.  2 

i 


tm 

*  S' 

z7T7~  ~r 

TTTV\ 

Cfc0SS-4*c.Ttoi4  eP 

UMMWUAt.  »P«ClMeM 

U -  w  - J 

fete  =  ~^C  fe< 


Describing  the  displacement  field  sketched  above,  we  assume 
it  to  be  of  the  following  form; 
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V  =  V(y) 


UJ  =  IAJ  (l)  (45) 


W 

At  x  =  /'I  there  exists  a  contour  for  u  which  describes  a 

curve  in  z 

U(!,4)  =  9  W  (46) 

2  2 

For  instance  g(z)  =  (z  -  h  )  yields  a  good  solution; 
however  for  now  we  take  a  more  general  approach  and  only 
assume  that  g(-h  )  =0. 

This  contour  is  continuously  modified  as  one  proceeds  to  the 
inner  regions.  This  will  be  represented  by 

U(*,i)=  (47) 

Since  the  displacement  in  the  coating  is  not  zero  at  the 
specimen  boundary  but  rather  produces  a  constant  strain  equal 
to  that  of  the  specimen,  another  term  must  be  added. 

=  400^(4)  4  C‘*  (48) 

The  other  displacements  are  chosen  to  give  constant  strains 


v(y)  -  C*y 

(49) 

10  il)  z  C-^i 

(50) 
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The  strains  become 


u-  ^ *  ci 

dx. 


s>  -  75  ■ c* 
<•  ••  ■  c> 


o*Y 

*X4 


2>U  .  ^vT  _  0 

^7  3? 

+  c>i*j  ,  -f  (*)  a\ i) 

"$*  " 

^VT  .  ^lO  -  O 

71  ^7  ’ 


The  stresses  become 


*  c,  »ct+c 


<?.,  -- 


— f  cj  4  ^U'q  t  C,  *  Ct*Cj] 


--  6  +  g' 


The  general  form  of  the  stress  expressions  is: 


<JT.  ^  P(K)6(i)  *  CoA^UAi 


Using  this  general  form  as  a  guide,  a  more  general  solution 
is  sought  using  the  field  equations.  Assume  TXy  *  'X2  =  0 
and  -s*  =0.  The  Equilibrium  Equations  simplify  to 

°  y 


'x  .  2>r, 

—  4* _ 


O  •  X4  „  O 

TT  ' 


4  5  O 

ST  5T 


These  equations  can  be  solved  with  a  stress  function  by 
prescribing  that 

(T  -  £l  =  -  tX  ^  -  ll  (55) 
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The  Compatibility  Equations  reduce  to 


W*  + 

7l(Ty  -- 
7lC1  + 
VlT  + 

V  Xt 


i  2>le 

l+i>t  T) ? 


rk  ^ 

_j_  7b 

l+0c 


o 


o 


o 


©  =  CTX  *  Gy  4-  (T* 


(56) 


Adding  the  first  three  compatibility  equations  gives 

f-*^)  Vz6  =  o 

(.  i  +>w 

Since  7l<Ty  =0 


Then 

l(<r*  +  <r.)  =  o 


(58) 


Substituting  the  stress  functions  into  (58)  ,  we  obtain  the 
governing  equation  for  it. 

V*  §  =  °  <: 

Using  the  biharmonic  theorem, 


§  ^  i<fi  +•  / 


(60) 


where  =  O  and  V L  -  O 

Using  separation  of  variables  for  ^  and  ^ 

^  =•  'f  (x)  ^  (i)  /■-  K(x)k(2) 

the  stress  function  can  be  represented  as 

£(x,i)  =  if(x)cj(£)  +■  k(x)k(2) 


(61) 


(62) 


The  stress  of  principal  concern  here  is  ^  . 

a;-  -foo  +  kOO  k"(i) 

d*1 

To  make  it  match  the  general  form  of  equation  (53) ,  let 

fw=  i  il  [19W]  --  c, 

<d  i* 


(63) 


(64) 


g(z)  is  found  by  simple  integration  and  h(x)  and  k(z)  are 
found  from  solving  uaPlace's  Equation: 


k(x}  =  A,  Co^K  Ax  +  A2  SuaK  Ax 
\cii):  6,  Cos  Ai  +  62  Sih  Ai 


(65) 

(66) 
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1 


where  ^  is  the  separation  parameter.  The  stress  function 
is  now  determined. 

5  Cx  i  +  C5  +  (A1CoiJ'\XX4  AtSinWXxJ^CosXl-j^  su  >*)(67) 

2. 

Differentiating,  the  stresses  become 

<JX  Z  C,  -  A^A,  Cos»K  Xx  +  A2  Si-nW  Cos  Xi  +■  ft*  S  i*.  Xt) 

r  A2 (A, Co^K  Xx  +  At  5«.aX>x)(6,  Co&  X*  +■  E>Su^  X2-)  (6  8) 

T  =  -Xl(AxSlwkXx  +  At  CoiK.>x)(7&,Suv>-t +•  ^CosXt) 


The  normal  stresses  must  be  symmetric  in  x;  therefore, 

A2  =  0 

Without  loss  of  generality,  A]_  =  1.  This  simplifies  the 
stress  expressions  somewhat  to  give, 

(jy.-C.  ~  A*C©skAx  {&,  Cos  Xi  +  S  L*\  Xi) 

A  I 

(j^  :  A*  C©sk  Ax  (£>,  Cos  Xi  +•  Sin.  A*) 

7:  =  Xu  SiyvkXx  (b,  Sun.  Ai  -  g^CosAi) 
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We  require  that  the  shear  stress  be  zero  on  the  free  surface 
z  =  0.  To  accommodate  this,!^  is  set  equal  to  zero.  The 
stresses  simplify  further  to  the  form 

<rK  -  c,  - 

( J  -  >"16,  C.Osk  CoS  Al  (70) 

i  t 

Going  back  to  the  compatibility  equations,  others  can  now  be 
satisfied.  To  start  with, 

Q  —  01  4  @y  4-  -  C-i  +  (71) 

Because  =  O  from  equation  (70),  equation  (56)  and 

(71)  give 

^  iiC'+^y)  =  =o  nil 

2>x.*-  o'*1-  7xv 

Integrating, 

(fy  =  *  ?U)  +  6(*)  (73) 
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Substituting  into  equation  (56) 


and 

(J  z  *  (C4  £  *  (£'+  1  V  Cl)  (77) 

Since  ^ ^ CT^  -O  from  equation  (70);  then  equation  (56) 
gives 

fe  -  i_  fc,  +  x(C<i*CS>  -  o  (78) 

Ti1 '  J 

This  is  satisfied  for  all  values  of  the  constants.  Since 
-  O  from  equation  (56) 

^  _  [c,  +■  )t(a€+Cs)  +'(^+'C7)  (79) 
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This  equation  is  satisfied  if  C4  =  0;  therefore,  all  of  the 
field  equations  have  been  satisfied  and 

(Ty  -  Cs* +  c,  iso) 

The  strains  in  the  plane  of  the  coating  are  now  of  the  form 

- 

(^z  -L  CobM- A:  (Cs  K+C+I  +  Cj) 

L 

(81) 

ts)  -  ~  f  Cff*  +  + 


The  loading  boundary  conditions  are:  ^-*.5  ;  ~  £ yj  A.t  £- 

This  gives 

(82) 

Mys^*~C^C7'^C‘l 


To  satisfy  these  requirements,  we  set 


A  z  JLl  cukxL  £$  z  O 

ZK 


(83) 
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pi 

* 


i 


Cg  is  set  to  zero  to  seek  a  solution  for  a  constant  cy. 

Since  -  —  Oj,  £y$  ,  equations  (82)  simplify  to 

-  \)b  £,c  -  C  |  -  C7 

(84) 

ec%s:  Ct  '  i><-c 1 

Solving  for  the  constants 

<*,  -  ec  tp  (4-->V) 

(85) 

C,  -  6-c  tp  (l  - 

l-Ti.” 


The  last  constant  is  found  by  setting  the  average  trans¬ 
verse  stress  equal  to  zero  at  the  boundary  of  the  specimen. 
Assuming  the  specimen  to  be  W  units  wide,  (J^so  ^  • 


From  which. 


aC  £c  tp 
i> -c.o*K.  (.?&) 


(87) 


The  stress  expressions  in  final  form  are: 


GI-- 

(j  .  (i-  OcOj)tyi 

6  -  ) 


(88) 


a  and  t  are  neglected  as  being  approximately  zero 
compared  to  As  can  be  seen  from  the  expression,  if 

vc  «  vg,  the  stress  state  becomes  uniaxial  in  the  limiting 
condition;  however,  if  there  is  a  Poisson's  Ratio  mismatch, 
and  there  normally  is,  a  transverse  stress  is  created  in  the 
coating  which  influences  to  some  degree  the  accuracy  of 
the  readings.  This  can  be  accounted  for  through  this  analy¬ 
sis  by  formulating  an  expression  for  an  equivalent  Poisson’s 
Ratio.  To  do  this  we  first  formulate  expressions  for  the 
average  strains  through  the  thickness  of  the  coating,  which 
will  be  what  the  fringes  will  be  measuring. 


i  - 


O-'W')  ] 


c  (iic) 
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An  equivalent  Poisson's  Ratio  is  defined  in  the  obvious  way. 


Substituting  elastic  material  properties  referred  to  earlier 


l.oo  ^4-  +■ 


.o  Shi 


.02.04 


C5£) 

c^k  (&) 

C-oik 


(91) 


A  plot  of  this  expression  was  made  assuming  that  W>>hc;  for 
instance ,  let  W  =  20hc>  A  plot  of  the  behavior  is  shown  in 
Fig.  7. 

If  the  specimen  material  is  taken  into  the  plastic  range  where 
'O.  c  0.500,  the  expression  for  Poisson’s  Ratio  becomes 


.  4-m  - 

.  1  Z 0~l 

- 

.  o4$S 

it) 

2_W<J 

Ccil'  (.T^) 


(92) 


For  the  same  parameters  of  W  *  20hc/  the  plot  of  equivalent 
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.3800  vt=  3793 


r 


0 

.C 


u 

-c 

fO 


o 

JC. 

CJ 


tO 


ro 


fO 


io 


fO 


.3200 


Poisson's  Ratio  for  the  plastic  case  is  given  in  Fig.  8. 

In  both  cases  the  combined  specimen  has  at  the  free  boundary 
the  Poisson's  Ratio  of  the  coating;  however,  this  falls  off 
rapidly  as  one  moves  to  the  interior.  Within  three  or  four 
coating  thicknesses  of  the  edge  the  Poisson's  Ratio  has 
approached  a  constant  value,  but  it  is  different  than  the 
specimen  value  alone. 

EXPERIMENTAL  COMPARISON  OF  STRESS  CONCENTRATION 
IN  A  STRIP  WEAKENED  BY  A  CIRCULAR  HOLE 
To  test  the  validity  and  accuracy  of  utilizing  the 
photoelastic  method,  a  test  case  using  a  circular  hole  speci¬ 
men  of  finite  width  was  selected.  R.  C.  J.  Howland  (8) 
determined  the  stress  distribution  around  the  circular  hole 
and  with  the  known  theoretical  solution,  it  was  possible  to 
determine  the  accuracy  of  the  photoelastic  method. 

EXPERIMENTAL  PROCEDURE 

A  specimen  was  designed  to  accommodate  a  2.0  inch  diameter 
hole  in  the  center.  To  determine  a  reasonable  width  and 
length  which  would  allow  a  uniaxial  field  to  exist  between 
the  center  hole  and  the  end  grips,  a  1/8  scale  PS-lC  model 
was  tested  under  various  tensile  load  conditions.  It  was 
determined  that  a  minimum  of  a  five-to-one  ratio  for  length- 
to-width  must  exist  to  eliminate  interference  of  fringe 
patterns.  Fig.  9  shows  the  geometry  selected. 
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The  aluminum  specimen  was  placed  under  an  8000  pound 
load  in  the  OTS  machine,  and  compensator  readings  were  taken 
with  a  Model  032  Photolastic  Compensator  at  the  90°  position 
on  each  side  of  the  hole.  They  were  found  to  be  122  /  121 
counts,  or  an  average  of  121.5  counts.  The  strain-optic 
equation  from  Appendix  A  can  be  written  as 

£y  “  ti  -  M  (93) 

At  the  circular  boundary,  -  -l)c 


4  -  P<r  Is*  _  o.  cob 2.9S  U/tk 

ty-  _  - 


Making  the  correction  for  reinforcement, 

£K-  \.OI  ly  -  O.OO-bbA* 


(94) 


(95) 


The  local  hoop  stress  is 


c r;':  w. 


Howland  * s  results  for  this  case  gave  34,719  psi  at  the  edge 
of  the  hole.  This  indicates  a  difference  of  2.4%.  Later 
we  will  show  that  this  difference  is  reduced  by  more  than 
a  factor  of  three  by  extropolating  to  determine  the  surface 


48 


stress  from  readings  in  the  interior. 

OBLIQUE  INCIDENCE 

In  order  to  determine  stress  values  away  from  the  edge 
of  the  hole,  oblique  incidence  was  used.  Fig.  10  is  a 
schematic  showing  the  coordinate  system  for  rotation  about 
the  Y  axis.  The  governing  equations  become, 


Figure  10.  Oblique  Incidence  Coordinates 
for  Rotation  about  the  y  axis . 


The  strain  in  the  x*  direction  is  given  from  the  strain 
tr ans  f ormat i on  laws 


i.%.  =  £*  Co>Le  +  6* 


(99) 


For  plane  stress,  the  following  stress-strain  relationships 
hold: 


t 


(100) 


(101) 


(102) 


Putting  these  together, 

-  .if-O  +  t>) 

l-Vc 


(10  3) 


and  the  expression  of  strain  along  the  oblique  axis  becomes 

t!  -  £,  to^Q-  ^  (U* 

K  t  X  <104) 

\  -V{, 
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From  Fig.  10  it  can  be  seen  that  Vvc  —  K.  (Let  &  .  Putting 
all  of  these  expressions  into  Equation  (98) ,  an  equation  for 
N  in  terms  of  known  quantities  can  be  obtained. 

- ' - ]Uv(l--i)cC«4t©)+**('X~C!>st£)]  1105) 

Solving  this  equation  simultaneously  with  equation  (97), 

*  r  [l  —  -  -  Nl(l-l)c  Co^©)|  (1G6) 

(ua)s-0  L  ;  J 

C  __  Nj©x)  (l-0c)C©*£  '  ^  ^  ~^c)  (107) 

^‘(1+Oc)s.^L  ]  J 


For  rotation  about  the  other  axis, 


frx  - 


£  (*/-*•) 


(108) 


Fig.  11  shows  a  schematic  of  this  oblique  incidence 
arrangement.  The  strain  in  the  oblique  direction  turns  out 
to  be  of  the  same  form  as  before. 


CoSe  +• 


(109) 
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Figure  11.  Oblique  Incidence  Coordinates  for 
Rotation  about  the  x  axis. 


Using  the  expresion  for  6.4  previously  derived  and  the 
expression  for  H  ,  the  expression  for  the  oblique  fringe 
value  becomes 


-  £b 
h 


M  = 
,Nex.  f 


6-0t)  Coi6 


Solving  this  simultaneously  with  the  normal  incidence 
equation  (97),  the  strain  expressions  are: 


Ft 


dK= - r 


N/’(Cosl©-Oc)  -  Me*(i-00  Cos©J 


die 


di] 


t\l  - 


Pt. 


M(i-v>cCc4v©)  -  ^,(i-i,.)Co4e] 


Calibrating  the  Polariscope 

In  oblique  incidence  the  determination  of  the  angle 
can  be  tedious  because  of  changes  in  it  at  the  surface  due 
to  the  index  of  refraction  being  different  in  the  coating 
than  it  is  in  air.  The  light  ray  is  bent  as  shown  in 
Fig.  12.  It  is  more  reliable  to  obtain  the  value  experimen¬ 
tally.  For  instance,  equations  (  106),  (107)  (111)  and  (112) 

can  be  written  as 

W-  Ft [am*, -cm] 

-  Ft  [aM*,  -  fcN] 
ex.  f *[e»M- 

1 

£y  Ft[cN  -  A 

Subtracting  (114)  from  (113)  or  (115)  from  (116), 

ty-U-  f^Ft(c-e>) 

Since  C  -  B  *  1  from  the  basic  equation  (A-l) , 

C  -  B  +  1 


(112) 


(113) 

(114) 

(115) 

(116) 

(117) 
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Figure  12.  Definition  of  Angle  of  Refraction 
for  the  Polariscope 
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Solving  two  of  the  other  equations  simultaneously  for  A  and  B 


A  -  L!2  - r  (118) 

(t y  -  6*)00«y 
£v  +  6k.  N^y 

6-  - - - — - - r  (ii9) 


Using  a  uniaxial  test  to  make  the  calibration  tests,  the 
expression  for  A,  B  and  C  become 


A=  O-o)  ^ _ 

( l  +  o)  (M^y  -  *0  ©*.) 


(120) 


isle*  -  3  N»y 
(t+o)(Mey- 


(121) 


c  -  *>  +  I 


(122) 


Several  readings  were  made  of  normal  incidence  and  both 

oblique  incidence  fringes  on  a  uniaxial  specimen  of  coating 

material  at  different  load  levels.  N.  and  N~  were 

ox  u  y 

plotted  versus  N  and  a  linear  regression  was  made  of  the  data 

=  0-S\0  f  O.  (123) 

-  0.  9969 
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N.,  -  '2166  +  1164  M  (12  4) 

fc.  -  0.9997 


Substituting  into  equations  (120)  to  (122)  and  using 
-0  =  9c.=  0.3793,  the  constants  were  found  to  be 

A  =  1.372 
B  =  0.899 
C  =  1.899 


Generalizing  for  Any  Material 

In  order  to  use  the  equations  of  oblique  incidence  for 
measurements  on  materials  of  any  Poisson's  Ratio,  the  ex¬ 
pressions  for  the  constants  were  used  to  solve  for  the  overall 
incidence  angle  ©. 

A-"0c)  Coi>G  _  I  t-jT. 

A  -  1.  J _ -  1  (125) 

(  ,  4-00 

Solving  this  equation,  0  =  31.85°.  The  other  two  expressions 
for  B  and  C  both  gave  0  =  31.76°.  The  good  agreement  gave 
confidence  in  the  numbers,  and  the  average  value  of  31.8° 
will  be  used  in  all  further  calculations. 

The  expressions  for  strain  in  terms  of  N  and  Ng^  become 
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£*  z 


Vi 


^=r  -(3.4,012.  -2.4,012. 0)mJ 

El.  2,.o«,oi  (l -a))  fsUy  -(z.4 ,oiz,  - 12.  0) Ni  ] 

1+0  L  -* 


(126) 


(127) 


For  rotations  about  the  x  axis, 


_  Pt_j^,.fc0,L-»,/.clto)N  -  3.0t..7(l-5)(vJ*,j  (128) 

"  t+o[;  J 


PHOTOELASTIC  DATA  FOR  THE  HOLE 

Normal  and  oblique  compensator  readings  were  taken  at 
0.0625  inch  intervals  at  the  90°  position  on  each  side  of 
the  hole  as  shown  in  Fig.  9  for  a  loading  of  12,000  lbs., 
or  a  far-field  stress  of  16,124  psi.  Six  different  compen¬ 
sator  readings  were  made  by  two  observers;  i.e.  each  observer 
made  three  readings,  and  these  are  listed  in  Table  VI.  The 
average  and  standard  deviation  are  also  listed. 
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TABLE  VI 


CIRCULAR  HOLE  COMPENSATOR  READINGS 
NORMAL  INCIDENCE 


X 

0.0625 

148 

151 

149 

149 

154 

14  8 

AVG. 

149.8 

STD  DEV 

2.32 

0.1250 

125 

130 

126 

124 

130 

125 

126. 7 

2.67 

0.1875 

108 

110 

108 

108 

108 

107 

108.2 

.98 

0.2500 

94 

96 

97 

98 

97 

97 

96.5 

1.38 

OBLIQUE  INCIDENCE 


X 

0.0625 

126 

125 

124 

125 

128 

124 

AVG. 

125. 3 

STD  DEV 

1.51 

0.1250 

102 

102 

100 

102 

104 

101 

101.8 

1.34 

0.1875 

85 

85 

84 

85 

87 

86 

85.3 

1.03 

0.2500 

74 

75 

75 

75 

76 

75 

75.0 

.63 

By  Chauvenet's  Criterion,  data  will  be  excluded  where 
probability  of  occurrence  would  be  less  than  ^  .  For  six 
readings  this  is  .0833.  Assuming  a  Gaussian  distribution 
of  errors  on  the  readings,  this  represents  1.73a  away  from 
the  mean  before  the  datum  becomes  suspect.  Table  VII  shows 
the  upper  and  lower  limits  of  acceptable  data. 
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TABLE  VII 


LIMITS  OF  ACCEPTABLE  DATA 


X 

NORMAL  INCIDENCE 

UPPER  BOUND  LOWER  BOUND 

OBLIQUE  INCIDENCE 

UPPER  BOUND  LOWER  BOUND 

0.0625 

153.8 

145.8 

127.9 

122.7 

0.1250 

131.3 

122.0 

104.1 

99.5 

0.1875 

109.9 

106.5 

87. 1 

83.5 

0.2500 

98.9 

94.1 

76.1 

73.9 

Data  falling  outside  of  these  ranges  were  eliminated  and 

new  averages  computed.  These  results  are  listed  in  Table  VIII. 


TABLE  VIII 

SMOOTHED  CIRCULAR  HOLE  COMPENSATOR  DATA 


DATA  POINT 

COMPENSATOR 

READINGS 

(in. ) 

Cn 

C9X 

0.06  25 

149.0 

124.8 

0.125 

126.7 

101.8 

0.1875 

107.8 

85.3 

0.25  0 

97.0 

75.0 

Extrapolating  to  the  Edge  of  the  Hole 

Because  it  is  difficult  to  measure  the  fringe  count 
right  at  the  edge  of  the  hole,  these  data  were  plotted 
and  extrapolated  back  to  that  point.  It  was  found  that 
log  log  values  of  the  compensator  readings  plotted  linearly 
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f 


with  x.  The  best  fit  was  obtained  by  using  the  first  three 
points. 

=•  0.3390  -  0-Z&S3*  (130) 

with  E  =  -.9996.  Using  the  value  for  x  =  0, 


C9x  -  152*31 

For  this  Compensator,  47  counts  constituted  1  fringe. 

89X  *  3'24 


Similarly  for  the  normal  incidence  fringes, 

Cn)  -  <5  «n  -  o.nzi* 


(131) 


with  R  =  .99996.  Using  the  value  for  x  =  0, 


and 


CN  =  176.81 

N  *  3.76 


At  the  edge  of  the  hole  the  radial  stress,  a ,  is  zero, 
which  requires  for  the  two  dimensional  plane  stress  geometry 
that 

lr  --  -  O  tQ  (132) 


Substituting  equations  (128)  and  (129)  into  equation  (132), 
a  relationship  is  found  between  and  N  for  all  “0  at 

8  X 
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the  hole  boundary. 


N**  =  0.  84  99  Nl  <133) 

After  extrapolation  the  measured  values  of  N  =  3.76  and 

Na „  =  3.24  check  out  to  within  1.3%  of  the  theoretical  values 
a  x 

from  equation  (133)  . 

Accounting  for  Poisson's  Ratio  Mismatch 

These  fringe  values  are  extrapolated  from  interior 
regions  where  "0  -  .3256;  however,  at  the  edge  of  the  hole, 

*0  =  .3793;  therefore,  they  must  be  adjusted  to  the  appro¬ 

priate  Poisson's  Ratio  value.  This  is  done  by  knowing  that 
at  the  boundary  the  hoop  stress  is  given  by 


Et  Ft  M 


l  +1> 


(134) 


For  a  given  ct0  existing  at  the  edge  of  the  hole,  the 
relation  between  the  photoelastic  readings  for  a  different 
Poisson's  Ratio  is  given  from  (134)  to  be 


_  Mj. 

I  +3,  l  + 


(135) 


This,  of  course,  assumes  E„  and  F  are  held  constant.  Thus 

s  _  e 

for  our  case,  0,  *  .3793  and  .3256,  and 

i 

IS|  1  2  i .  040$  N 
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The  actual  normal  incidence  fringe  value  at  the  edge  of 


the  hole  is  now 


M  -  (l.O40s)(W?)  =  3.91 


To  find  the  corresponding  value  for  the  oblique 
incidence  fringe,  equation  (133)  is  used  in  conjunction 
with  equation  (135) . 


r  1.D40S  Ne* 


»=  .vm 


(136) 


From  this  the  value  of  Nfl  at  the  edge  of  the  hole  is 


Ne*  =  (i-o*os)(i.u)=  5*7 


The  last  correction  to  the  data  is  for  reinforcement, 
which  for  these  materials  is  a  factor  of  1.0162  from 
equation  (39) . 


N  =  3.98 


NQx  =  3.43 


Check  against  Numerical  Solution 


A  check  can  now  be  made  against  the  numerical  solution, 

which  for  this  geometry  and  loading  is  given  by  a  finite 

element  solution  developed  by  Kaiser  (4)  to  be  c.  =  51,678  psi. 

0 

Using  N  *  3.98,  the  extrapolated  photoelastic  measurements 
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give  from  equation  (134). 


^  .  si.g“P»w 

®  1.  *.193 

This  value  is  within  0.3%  of  the  finite  element  value. 

This  is  surprisingly  good  agreement. 

Check  of  Stress  Distribution 

A  check  of  the  accuracy  of  the  photoelastic  work  was 
also  made  on  the  stress  distribution  away  from  the  hole  by 
means  of  a  finite  element  solution.  The  criteria  established 
to  determine  uniform  boundary  stress  distribution  was  uni¬ 
formity  in  nodal  displacements  along  the  loaded  edge  as 
discussed  by  Segerlind  (9) .  In  the  models  used  for  FINITE 
ELEMENT  ANALYSIS  (FEA)  in  this  report,  nodal  displacements 
were  uniform  to  within  0.1%,  and  the  resulting  stress  distri¬ 
bution  was  uniform  axially  to  within  0.1%  at  the  panel  ends. 

In  all  cases  two-dimensional,  eight  nodea,  isoparametric 
elements  were  used.  These  higher  order  elements  cannot  be 
loaded  in  an  "intuitive"  manner,  as  discussed  by  Zienkiewicz 
(10).  Fig.  13  shows  the  nodal  loading  required  to  obtain  a 
uniform  surface  load. 

Element  Meshes 

Two  meshes  were  developed  for  each  panel  analyzed. 

A  reasonable  effort  was  made  to  keep  element  corner  angles 
as  close  to  90°  as  possible  to  reduce  the  effect  of  element 
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distortion  discussed  by  Hopkins  and  Gifford  (11) .  All 
meshes  modeled  a  quarter  of  the  actual  panel  by  using  the 
two  axes  of  symmetry  as  is  common  practice  in  FEA.  The 
step  from  course  to  fine  element  meshes  was  made  so  that 
each  element  in  the  course  mesh  was  subdivided  into  four 
smaller  elements  of  the  same  type.  Such  a  mesh  subdivision 
can  be  expected  to  give  monotonic  convergence  of  results, 
as  shown  by  Cook  (12) ,  and  allow  extrapolation  to  results 
of  an  infinitely  fine  mesh.  Figures  14  and  15  illustrate 
the  element  meshes  used  in  this  analysis. 

Using  ADINA 

Once  the  mesh  was  developed,  input  data  was  prepared 
in  accordance  with  the  ADINA  user's  manual  (12) .  This  same 
set  of  data  was  used  as  input  for  PSAPl,  the  plot  routine,  to 
check  for  errors  and  provide  a  graphical  display  of  the 
element  mesh.  After  preprocessing  by  PSAPl,  the  data  is 
entered  into  ADINA  for  analysis.  In  the  case  of  linear 
analysis,  two  types  of  stress  output  may  be  specified, 
nodal  point  or  Gauss  integration  point.  Nodal  point  output 
can  be  computed  for  up  to  eight  node  point  stresses  for 
each  element.  Since  2x2  Gauss  integration  was  used,  four 
Gauss  point  stresses  were  computed  for  each  element.  The 
2x2  Gauss  integration  is  recognized  as  the  most  efficient 
integration  order  for  this  type  of  analysis  [Ref.  10,  p.  284], 
The  linear  analysis  used  an  isotropic  linear  elastic  material 
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FIGURE  15 


FINE  MESH  FOR  CIRCULAR  HOLES 
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model  (MODEL  "1"  in  Ref.  4)  which  required  input  of  tensile 
modulus  and  material  properties. 

For  static  analyses  ADINA  uses  a  time  function  method 
to  apply  loads  in  steps.  Linear  analysis  loading 
was  accomplished  in  a  single  step  to  a  nominal  value  of 
3,000  lbs  load.  The  stress  output  from  ADINA  is  a  listing 
of  nodal  or  Gauss  point  stresses  for  each  element.  Since 
the  only  area  of  interest  in  this  analysis  was  the  distri¬ 
bution  of  stresses  along  the  reduced  cross-section,  no  large 
post-processing  program  was  developed  or  used.  All  final 
computations  using  ADINA  output  data  were  accomplished  on 
a  HEWLETT-PACKARD  9 8 30 A  calculator,  using  short  programs 
coded  in  BASIC.  If  more  extensive  stress  distribution 
information  were  desired,  some  form  of  automated  post¬ 
processing  would  be  necessary  to  reduce  the  computational 
workload.  At  a  minimum,  nodal  stress  outputs  by  ADINA  must 
be  averaged  to  obtain  unique  values  of  stress  at  nodes  shared 
by  more  than  one  element. 

Richardson  Extrapolation 

The  use  of  course  and  fine  meshes  allows  extrapola¬ 
tion  to  an  infinitely  fine  mesh  by  means  of  the  Richardson 


extrapolation  (14)  where: 

tfc(Kp)  -  CTf(Kc) 


Gjuirx p 


kp  -  K 


(137) 
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where 


r 


sional  analysis.  Better  extrapolations  can  be  obtained  by 
using  three  or  more  defined  meshes,  but  the  computational 
effort  increases  significantly. 

Optimal  Stress  Locations  and  Local  Smoothing 

It  is  generally  accepted  that  the  most  accurate 
sampling  points  for  stresses  are  the  Gauss  integration 
points  within  the  element  [Ref.  10,  p.  281].  In  this  analy¬ 
sis,  the  nodal  points  are  of  the  greatest  interest;  thus 
a  technique  of  local  smoothing  must  be  applied  to  the  inte¬ 
gration  point  stresses  to  obtain  nodal  stresses  as  reported 
by  Hinton  and  Campbell  (15) .  The  nodal  values  obtained 
must  then  be  averaged  if  shared  by  two  or  more  elements. 
Circular  Holes  in  Linear  Material 

The  FEA  results  for  a  circular  hole  in  a  finite  width 
strip  were  used  to  validate  the  photoelastic  procedures  dis¬ 
cussed  earlier.  The  results  of  Howland  (8)  were  compared  to 
both  the  Gauss  point  smoothed  results  and  the  nodal  output 
results  in  Fig.  16.  The  stress  concentration  factor,  o/am, 
is  referenced  to  the  far-field  stress.  The  smoothed  results 
give  the  best  match  to  the  results  of  Howland  at  the  edge 
of  the  hole,  and  the  only  significant  variation  between  the 
two  FEA  methods  occurs  within  the  first  0.25  inches  from  the 
edge.  In  order  to  obtain  the  0.25  inch  stress  value  for  the 
coarse  mesh,  in  the  Gauss  point  smoothed  result,  a  midside 
node  value  had  to  be  obtained  by  the  averaging  method. 
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FIGURE  16 

CIRCULRR  HOLE  \^.2  LINEAR  RESULT? 


DISTANCE  FROM  HOLE  /  inches 


TABLE  IX 


X  =  0.25  FEA  RESULTS  -  NODAL  OUTPUT 


DISTANCE 

PROM  HOLE,  in. 

a oo  «  4500  psi 
<J0  *psi  ar  ,psi 

0.00 

14745.5 

329.8 

0.25 

9447.0 

2181.1 

0.50 

7232.6 

1599.3 

0.75 

6304.7 

1473.5 

1.00 

5771.7 

1139.6 

1.25 

5471.2 

891.6 

1.50 

5266.9 

671.0 

1.75 

5091.7 

489.1 

2.00 

4940.2 

354.5 

TABLE  X. 

DISTANCE 

FROM  HOLE,  in. 

X  =  0.25  FEA  RESULTS  -  GAUSS  OUTPUT 

a,,,  =  4500  psi 
a0  ,psi  ar  ,psi 

0.00 

14422.8 

121.2 

0.25 

8721.2 

1927.1 

0.50 

7317.1 

1639.7 

0.75 

6218.6 

1452.4 

1.00 

5779.6 

1151.5 

1.25 

5453.3 

891.0 

1.50 

5267.1 

677.3 

1.75 

5085.9 

486.9 

2.00 

4937.9 

343.8 
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The  linear  distribution  of  smoothed  stresses  along  the 
sides  of  the  element,  (15) ,  appears  to  produce  a  less  accurate 
result  in  this  area  of  extreme  stress  gradient,  when  compare! 
to  ADINA's  nodal  output  result.  This  tendency  was  noted 
in  all  cases;  however,  the  peak  stress  values  from  smoothed 
results  consistently  gave  better  correlation  with  other  inves¬ 
tigators  (16)  . 

A  circular  hole  with  X=0.25  was  also  analyzed  and 
compared  to  the  photoelastic  experimental  results.  The 
stress  distribution  is  given  in  Table  IX  and  Table  X.  The 
data  for  the  first  two  locations  has  been  normalized  and 
entered  into  Table  XI  to  show  a  comparison  with  photoelastic 
results  similarly  normalized. 


TABLE  XI 

COMPARISON  OF  STRESS  DISTRIBUTION 

Finite  Element  Photoelastic 


X 

C  j/lJ« 

Q 

\ 

Q 

8 

ae/c. 

0.00 

-.0269 

3.2051 

-.1167 

3.0780 

0.25 

.4282 

2.0993 

.4111 

2.1627 

Comparisons  of  the  photoelastic  data  with  the  finite 
element  show  agreement  within  3%  and  4%  on  the  average.  This 
is  not  as  good  as  the  correlation  made  at  the  edge  of  the 
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hole  where  we  required  or  =  0.  If  c  is  allowed  to  be 
whatever  the  readings  show;  then  the  correlation  is  in  the 
3%  or  4%  range;  however  as  can  be  seen  at  x  =  0 ,  the  finite 
element  is  not  exact  either,  especially  at  these  small  x 
values . 

Sensitivity  Studies 

At  the  edge  of  the  hole  where  1?  =  .3793,  the  expres¬ 
sions  for  strain  from  equations  (128)  and  (129)  become 


-  I.  S198  xi©"'**  ”  Z.AL9io  * 


(140) 


u  z  3.M*6*»* 


-  Z.42.%*  i©’*  ^ 


(141) 


Substituting  these  expressions  into  the  plane  stress  consti¬ 
tutive  equations 


(j\  z  Z.08^1  xio4M  -  Z .  4531  x  i  o4,  /vie* 


(142) 


0*  z.  Z-& ~~  Z.4S& 7xio 


(143) 


The  sensitivity  of  the  stress  is  given  by 

2W 


(144) 
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Applying  this  to  equations  (142)  and  (143) 


(i<rr  -  Z*089l  -  Z.4S&1  **04clhl9-K  (145) 

s  Z  &9$o  *  lO^/sI  -  Z.4-S&1  *s  10 4  <±tieiL  (146) 


Selecting  the  maximum  deviation  of  N  and  N  as  1.73a,  where 

0  X 

a  is  the  standard  deviation  of  the  smoothed  data 


dM  =  t.22)  r  ±2-11 

<4nU=  (i.7J.)(t  .m)  = 


Then  da  can  have  two  values  depending  upon  how  the  plus  and 
minus  signs  couple  up.  The  worst  of  these  are: 

dtrr  =  i  u.97 

d<re  :  ±  P5"' 

These  figures  represent  the  maximum  errors  that  would  be 
expected  due  to  misreading  the  fringe  pattern,  and  they  are 
approximately  4%  on  o9  and  explains  how  non- zero  results 
for  cr  arise. 


CONCLUSIONS  AND  PFSULTS 

From  very  careful  experimental  procedures  where  no 
handbook  values  for  material  properties  were  used  but  all 
were  carefully  determined,  photoelastic  measurements  checked 
theoretical  results  within  4%.  This  error  included  all  errors 
in  material  properties,  in  operator  readings  and  in  equip¬ 
ment  calibration.  Without  very  careful  and  detailed  measure¬ 
ments  on  every  parameter,  photoelastic  results  could  easily 
be  in  error  by  15%  to  20%. 

With  this  level  of  accuracy  established  by  controlled 
laboratory  procedures  and  calibrations  of  equipment  in  place, 
we  will  now  proceed  to  study  the  elastic  and  plastic  behavior 
of  notches  where  fatigue  failure  originate. 
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APPENDIX  A 


Equations  of  Photoelast:  city  Coatings 


( A-l) 


(A-2) 


(A- 3) 


(A- 4) 


( A-5) 
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